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Several authors have argued that self-consistent f{R) gravity models distinct from ACDM are 
almost ruled out. Confronting such claims, we present a particular two-parameter f{R) model that: 
(a) is cosmologically viable and distinguishable from ACDM; (b) is compatible with the existence 
of relativistic stars; (c) is free of singularities of the Ricci scalar during the cosmological evolution 
and (d) allows the addition of high curvature corrections that could be relevant for inflation. 
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Introduction. Since the discovery of cosmic acceler- 
ation, more than a decade ago, considerable effort has 
been devoted in cosmology to understand what is the 
physical mechanism responsible for it. A relic cosmologi- 
cal constant A, even though arguably the simplest expla- 
nation and in good accordance with observations, faces 
some theoretical difficulties (mainly due to the cosmic co- 
incidence problem and related fine-tuning [H) that have 
motivated an intense search for alternatives. These can 
be divided into two main conceptual approaches, both 
involving the introduction of new degrees of freedom (see 
for instance either one modifies the left hand side of 
Einstein's equations (modified gravity) or one adds a new 
term to the energy momentum tensor, arguably associ- 
ated with a new fundamental field not directly related to 
gravity. 

Special attention to the former approach has been 
given in the last five years. In particular, /(i?) grav- 
ity theory, due to its sinipHcity, received the main focus 
(for a recent review, see [3|| and references therein). This 
approach amounts to writing the action as 
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f{R) + Cn 



(1) 



where f{R) ~ R+A{R), R is the Ricci scalar and A(i?) is 
an arbitrary function. General Relativity (GR) without 
a cosmological constant is obtained in the special case in 
which A(i?) is identically zero. Although a great deal 
of effort has been employed to develop this approach, 
it appeared to be a difficult challenge to build a new 
Lagrangian that does not spoil the successes of GR — 
one that passes solar system tests, describes the early 
universe, allows a matter-dominated phase followed by 
an accelerating attractor j3] — and, at the same time, do 
not suffer from curvature singularities The presence 
of singularities may have devastating consequences and 
could forbid, for instance, the formation of relativistic 
stelar objects such as neutron stars 0. 

Singularity- Free f(R) Model. Several popular f{R) 



models investigated in the literature are generalized by 
the following expression 



fiR) = R-Rsl3ll 
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For instance, choosing /3 = — 1 we obtain the models 
presented in [7|| ; for /3 = 1 we recover the model proposed 
in for n = 2 we get the f{R) function discussed in j^. 
In this letter we consider the special case in which n = 1 
and we take the limit /3 — > 00. In this limit ^ can be 
recast as (rewriting Rs as ai?*) 



/(i?) = R- aR^\n{l + — 



(3) 



where a and i?* are free positive parameters. Notice that 
the above function satisfies the stability conditions [lo| : 
(a) fRR <ff/dR^ > (no tachyons [Hi); (b) /i?j= 
d//di? > (no ghosts) for a < {R/R„ -\- 1), where R is 
the value of the Ricci scalar at the final accelerated fixed 
point; and (c) limfl^ooA/i? = and Ihnji^oo = 
(GR is recovered at early times). Above and henceforth, 
Ar := dA/di?. 

Starting from the action |[T]) , one obtains the equation 
of motion for f{R): 



OfR - \f ) .9^. = SttGT^,, (4) 



the trace of which is given by 



□/H = ^T+i(2/-/Hi(), 



• 

where T is the trace of the energy-momentum tensor. 
We now introduce the scalar degree of freedom (d.o.f.) 
X /i? and write the equations based on the mapping 
from ,f{R) gravity (with positive first and second deriva- 
tives) onto Brans-Dicke scalar-tensor theory with param- 
eter w = 0. The resulting field equation is 
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Figure 1: V{'X)IR-, :— —U{x)/R* for different models: ours, 
with a = 2 (blue solid line), Starobinski's for {n=2, A = 
2} [see [^] (red, dashed line) and Hu & Sawicki's for {n=2, 
= 1, ci/c2 = 2} [see 0] ( green dot-dashed line). The 
physically interesting region is < x < 1- For the multi- 
valued potentials only the lower lines are physical. 



with the force term given by T :~ — (SttG/S) T and 
dV{R{x)) 1 



(2/ - fRR) . 



(7) 



When applying the model ([3]) to a spatially homoge- 
neous and isotropic universe the scalar d.o.f. becomes 



xm)] 



R{t) + i?* 



(8) 



and the d'Alembertian in ^ is effectively just a time- 
derivative: □ = —d'^/dt^ — 3Hd/dt; our choice for the 
metric signature is (—,+,+,+). It is straightforward to 
see that x^l~ asi?— >cx), which points out the 
same singularity jH| featured in previous models (sl. lol. [l^ . 
Inverting the relation ([8| and integrating ([7| we find that 
(up to a constant) 



-a(2x-3) In 
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X-3 



(9) 

Note that since ^ defines a one-to-one relation between 
X and R, the potential V{x) is well-defined and not multi- 
valued, contrary to the models in 0,S,@|. Fig. [1] depicts 
the potential for a = 2, as well as typical potentials de- 
rived from models 0, Taking the limit x 1~ we 
find that 
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which shows the presence of an infinite barrier at x = 1 
that prevents the singularity discussed in to be 
reached. 

We can understand this result in a more intuitive way 
by making use of the well-known duality between f{R) 
and scalar-tensor theories: a conformal transformation of 



the metric can cast the Lagrangian from the Jordan into 
the Einstein frame, converting the scalar d.o.f. x into a 
canonical scalar field x '■= — ■\/3/167rGlnx The field 
equation for x has the same structure of ([6]), but with 
the following potential 



Ve (i?(x)) 
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Afl)2- 



(11) 



All the discussion above, regarding the presence of an 
infinite barrier, applies to Ve as well. Note that since 
1 + Afl > (stability condition (b)), the numerator of 
Eq. (fTTj) is the only factor that can make the potential 
diverge as i? oo. In [lj| the singularity was avoided 
by introducing an extra high-curvature term ai?"'^^' in 
the model investigated in [9|. It is easy to see why that 
kind of correction works: in that case, the numerator 
in (fTTj) is itself proportional to i?". Nevertheless, such 
term cannot be used, at the same time, both to avoid 
the singularities and to generate infiation This is 

not the case of the model investigated in this letter, since 
it is not necessary to include such terms to avoid the two 
singularity problems, as we have shown above (for the 
case discussed in [H|) and will show bellow for the case 
discussed in [^. 

Notice that two different singularity-free classes of 
f{R) are possible: we can pick a function A such that 
either lim^^tx) RAr = oo or limfl_>oo A = — oo holds. In 
the former case, A can even become constant as i? ^ oo 
— which actually happens in the models previously men- 
tioned [1, 0| — but it should do so slowly, thus keeping 
the divergence of RAji, which does not happen on those 
models. The model |[3]) belongs to the latter case. An- 
other interesting example of this class is 



f{R) ^R-aRjl 



R 



(12) 



with a > , i?* > and n e (0, 1). Although preliminary 
tests indicate that this model is cosmologically viable, 
it carries an expHcit positive cosmological constant, in 
direct contrast to 

We further remark that the potential (flTI) derived 
from ^ generates a Yukawa-like force which is fully com- 
patible with the Chameleon mechanism [l5| . In other 
words, the mass of the x field is large (small) when the 
background matter density is large (small). This mass 
dependence on the local environment explains how this 
extra (or fifth) force can have cosmological implications 
while at the same time evading detection by local gravity 
experiments. 

Relativistic Stars. The authors of Q argue that the 
very existence of relativistic stars poses a strong con- 
straint on f{R) gravity theories. For the models studied 
in that paper, it was not possible to evolve the metric 
from inside a star up to large spatial scales and match the 
de Sitter solution asymptotically. We show below that 
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ds^ = -N{r) dt^ + 
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dn^ 
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Figure 2: The field x with a — 1.2, pc — O.Spo and Rc varying 
from 1 X 10~*po to 4 x 10~*po. The arrow points out the 
solution that stops at the maximum of the potential at r ^ 
oo. The (red) thin lines indicate the region inside the star. 



this divergence is circumvented by model ^ and, there- 
fore, does not represent a general feature of f{R) models. 
For the sake of clarity, we follow the classical-mechanics 
analogy used in and the necessary definitions. We con- 
sider a static and spherically symmetric metric and write 
the d'Alembertian in © as □ = d^/dr^ + {2/r)d/dr 
in spherical coordinates; we are assuming a Minkowski 
background for a moment. In this case, Eq. ([6]) can be 
seen as the equation of motion of a classical particle of 
unit mass (albeit one whose "time" coordinate is our spa- 
tial coordinate r) submitted to both an external and fric- 
tional forces. Therefore 



d^ 
dr^ 



2dx 
r dr 



(13) 



where J- := —J- and Tu := —AU/dx are, respectively, 
the "force" due to the trace of the energy- momentum ten- 
sor (non vanishing inside the star) and the "force" due 
to the potential U{x) = -^(x)i see Eq. ^ and Fig. [H 
Again, the change in sign is just a consequence of the fact 
that now it is the spatial (instead of time) dependence of 
X which is the most relevant. 

For the models analyzed in [^, there was no solution 
which would describe a particle going uphill pulled by 
the force T (while still inside the star) and stop at the 
top of the potential at r ^ oo, which would correspond 
to the de Sitter metric. The particle would either return 
and reach the singularity at x = 1 (where R — > oo) or 
overshoot the potential towards x = (which would also 
lead to a singularity, for instance, in the Kretschmann 
scalar K := R"^^" RaPnv)- Fairly enough, U{x) diverges 
at X = 1) as in all other models 0, @|. As we will show 
below, the advantage here is a well-behaved solution fully 
compatible with relativistic stars embedded in a de Sitter 
universe. 

Let us now determine the full evolution of the x field. 
As previously mentioned, we start from a static and 
spherically symmetric metric 



B{r) 

and assume a constant energy-density star 
whose energy-momentum tensor is given by 
= d\&g{—pQ,p{r),p{r),p{r)). The initial condi- 
tions at Ti = 10~®i?* i.e, close to the center of the 
star, are given by N{ri) = \ + N2rf, Biji) = \ + B2rf , 
p{ri) = Pc+V2 rf/2 and x{n) = Xc (l + C2 r}/2). 
The coefficients A^2, ^2, P2 and C2 can be written in 
terms of po = 2 x lO^Aoff and of the central values 
Pc = 0.3po, Rc = 10-Vo, V{xc) and AV/dx{Xc). 
The effective value of the cosmological constant is given 
by Aoff ~ -Ri/4, where i?i is the value of the Ricci 
scalar when AV/Ax ~ 0. We refer the reader to the 
original paper Q for the full set of equations. Energy 
conservation provides an important relation between 
p(r) and N{r) inside the star. We evolve the system 
{p, i?, X, dx/dr} from up to the radius TZ of the 
star (defined by p(7^) = 0) where we require continuity 
of the variables. From then on we evolve the system 

— 1/2 

{iV, X, dx/dr} until r = i?* (cosmological scales). 

We show in Figure [2] the behavior of the field x for 
different values of initial conditions. Note that some tra- 
jectories do not get past the top of the potential and 
return towards the singularity at x = 1 (top 3 curves) 
while others (3 lowest ones) overshoot and go towards 
X = and one (indicated by an arrow) stops right at 
the maximum. We recall that this solution was obtained 
without any high-curvature correction. It is obviously an 
issue of fine tuning the initial conditions to stop exactly 
there. Another remarkable feature of this model is the 
absence of singularity in if as x decreases below the peak 
of its potential. 

A Promising Model. A viable cosmological model 
must start with a radiation-dominated universe and 
have a saddle point matter-dominated phase followed 
by an accelerated epoch as a final attractor. We can 
formally state such criteria if we use the parameters 
TO := Rf^RR/fR and t := -R/r/J ■ We refer the 
reader to the original paper \^ for a full discussion 
on this subject. An early matter-dominated epoch of 
the universe can be achieved if m{x w — 1) w O"*" and 
dTO/dr(r w —1) > —1. Furthermore, a necessary condi- 
tion for a given model to reach a late-time accelerated 
phase is < TO(r w —2) < 1. The model ^ satisfies 
both constraints for a > 1 regardless of i?* . 

Using Q, we obtain the modified Einstein's equations 
below for a homogeneous universe filled with matter en- 
ergy density pm (baryons and cold dark matter) and ra- 
diation energy density pr- 



= SttG (prn +Pr) + {/rR - /) /2 - 3HfR + 
+ 3H^{l-fR) 
-2H = SttG ip,n + 4p,./3) + fR - HfR - 2H{1 



(15) 
(16) 
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Figure 3: Cosmological evolution of the densities fi,„, ^r, 
(solid lines), the deceleration factor q (dot-dashed line), the 
jerk j (dotted line) , the equation of state parameters Wx and 
Wes (dashed and dotted lines, respectively), for a = 2. 



where a dot corresponds to derivative with respect to t, 
H = a/a and a{t) is the scale factor. From the equations 
above, we can define px, Px and Wx ■= Px/px, respectively 
the energy density, pressure and the equation-of-state pa- 
rameter of the so-called "curvature fluid": 

SnGpx ifRR -f)/2- 3Hfn + 3H\1 - Jr) (17) 

87:Gpx fn + 2Hfn - {2H + 3H^)il - fn) + 

+ (/ - fRR)/2 . (18) 

These definitions are such as to guarantee that the cur- 
vature fluid is conserved and only minimally coupled to 
matter and radiation [l3|- We also deflne the relative 
densities fi^ (where i stands for either radiation, matter 
or curvature) fli := SnGpi/^H^. 

In Figure [3] we plot the behavior offlm, ^r, ^x, the de- 
celeration parameter q := —da/a^, the jerk j aa^/a^, 
the equation of state parameters for the curvature fluid 
Wx and for the effective fluid WeS := ptot/ptot = (pr + 
Px)/{Pm + Pr + Px) , all of which can be written in terms 
of known variables R, and pi. In Fig. [3] we can clearly 
distinguish the radiation-dominated era when g ~ 1 (and 
j ~ 3, not shown), followed by a transient domination 
by matter {q ~ 1/2 and j ~ 1), the current acceler- 
ated expansion {q < 0) and the final de Sitter attractor 
{q = —j = —1). We find similar results for different ini- 
tial conditions and parameters, indicating what seems to 
be an absence of fine tuning. We remark that the Wx 
curve in Fig. [3] is noisy in the early universe since at that 
time Px is too small and the numerical calculation of Wx 
becomes inaccurate. 

We point out that there is some residual arbitrariness 
in defining px and px even if one is only interested in con- 
served and minimally coupled fiuids. The one we follow, 
together with the definition of fij , is convenient for com- 
parison with GR-based interpretations of observations. 



As another consequence of Eqs. I|17p and lfT8|). Wx nei- 
ther crosses —1 nor diverges at low redshift in contrast 
with [l3|, where slightly different definitions are adopted. 
Note, however, that observable quantities like H and Pm 
are well-defined and in fact, using either definition, have 
the same cosmological evolution. 

Conclusions. We have shown that some recent results 
in the literature regarding divergences in f{R) theories 
are not as general as previously thought. In fact, even a 
compact, two-parameter lagrangian like the one in ([3]) 
can evade the aforementioned singularities. Observa- 
tional constraints on this model are under investigation 
and the results will be published elsewhere. 
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